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. Abstract This is mainly a brief review of some key achievements in a "hot" area 

^ . of theoretical and mathematical physics. The principal aim is to outline the basic 

structures underlying integrable quantum field theory models with infinite- dimensional 
symmetry groups which display a radically new type of quantum group symmetries. 
Certain particular aspects are elaborated upon with some detail: integrable systems 
of Kadomtsev-Petviashvili type and their reductions appearing in matrix models of 
strings; Hamiltonian approach to Lie-Poisson symmetries; quantum field theory ap- 
proach to two-dimensional relativistic integrable models with dynamically broken con- 
QhI formal invariance. All field-theoretic models in question are of primary relevance to 

diverse branches of physics ranging from nonlinear hydrodynamics to string theory of 
fundamental particle interactions at ultra-high energies. 



^ ! 1. Introduction 



One of the prevailing views in modern theoretical physics is that fundamental laws of 
Nature can be derived and understood in terms of field-theoretic models in a lower 
dimensional space-time possessing infinite- dimensional symmetry groups and, thus, as 
a rule being integrable. 

These models in their various facets and disguises are encompassed in the extremely 
rich and rapidly developing branch of string theory It is widely believed that string 
theory is the most viable candidate for a unified theory of all fundamental interactions 
at ultra-short distances which, in particular, will provide a consistent reconciliation 
between General Relativity and Quantum Mechanics - one of the major challenges of 
this century's Physics. 

We have in mind two large classes of integrable models : conformal field theory 
(CFT) 0, ^ and massive completely integrable models 0, H, || in = 2 space-time 

"'^To appear in "Mathematical Physics Towards the XXIst Century", R.N. Sen and A. Gersten eds. 
(Proc. of the Int. Conf. held at Ben-Gurion Univ. of the Negev, Beer Sheva, March 14-19, 1993) 

^On leave from the Institute for Nuclear Research and Nuclear Energy of the Bulgarian Academy 
of Sciences, Tsarigradsko Chausee 72, BG-1784 Sofia, Bulgaria; E-mail: dstoyan@bgearn.bitnet; Fax: 
359-2-755019 . 
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dimensions. Typical examples of CFT are the rational CFT's, of which the most ex- 
tensively studied are the Wess-Zumino-Novikov-Witten (WZNW) models for various 
Lie groups G and models, obtained from them by gauging of different subgroups H of 
G 0. Thoroughly studied models in the second class are : Sine-Gordon, (nonabelian) 
massive Thirring models, Toda models for various groups, the Korteveg-de Vries (KdV) 
and Kadomtsev- Petviashvili (KP) integrable soliton evolution equations and their hi- 
erarchies. The feature of integrability, common for both classes of models, stems from 
the infinite-dimensional Lie-algebraic structure they share : in the first class this being 
the Noether symmetry algebra, and in the second, the Hamiltonian structure. The un- 
derlying infinite-dimensional symmetries, manifesting themselves through the Virasoro 
(conformal) and affine Kac-Moody algebras [Q, as well as through various (infinite- 
dimensional) generalizations thereof - e.g. the VT-algebras 0, play a crucial role in 
solving and interpreting integrable models. 

The interrelation between CFT and completely integrable models became recently 
explicit through the appearance of KdV and KP integrable hierarchies in the matrix 
model description |]lOl of (sub) critical strings (i.e., 2-dimensional gravity interacting 



with conformal "matter" fields). Thus, it is precisely the integrable field theories which 
provide the proper framework for incorporation of the huge symmetries of string theory 
models. 

The property of integrability is studied in two aspects : (1) classical and (2) quan- 
tum. The principal questions, which one should answer, are the following : (la) 
classification; (lb) action- angle type of variables; (Ic) exact integration of the equa- 
tions of motion — in the classical case, and (2a) classification; (2b) exact scattering 
amplitudes; (2c) exact correlation functions of local fields off the mass shell — in the 
quantum case. 

The concepts and tools to approach the above topics invoke, besides the theory of 
infinite-dimensional Lie algebras and groups, another outstanding branch of mathemat- 
ical physics - symplectic geometry or, equivalently, Hamiltonian mechanics. The aim 
here is to uncover the geometrical foundations of the relevant field theories. The use of 
geometric techniques offers powerful means for unifying various physical theories and ob- 
taining new insights. It is just in the realm of string theory and the theory of completely 
integrable systems where the intertwining of Hamiltonian and Lie-group structures in 
field-theoretic models attains new immensely important qualities. Their quintessence is 
manifested in the development of the principal methods to solve the quantization prob- 
lems in integrable models : quantum inverse scattering method [|11], 0, representation 



theory of infinite-dimensional Lie algebras , Quantum Groups (non-commutative and 



non- CO commutative Hopf algebras) p!3|, [12 



The generic integrable models are massive field theories which in a sense can be 
regarded as integrable perturbations of conformal field theories |]T1|. Such models have 
the advantage of being relativistic invariant and classifiable by the conformal models 
of which they are perturbations. The latter describe the renormalization group fixed 
points of these massive integrable models. Their most essential feature, explicitly expos- 
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ing the intimate connection to conformal models, is the existence of multi-Hamiltonian 
structures, i.e. the existence of at least a second Hamiltonian structure which is com- 
patible with the canonical i?- matrix Kirillov-Kostant structure The corresponding 
fundamental Poisson brackets (linear i?-matrix brackets and quadratic (Sklyanin) R- 
matrix brackets) naturally arise and are exhaustively understood within the classical 
"inverse scattering" method 0. Also, they can be deduced in the semiclassical limit 
from the basic algebraic structures : 

(1) Yang-Baxter equation for the quantum version of the i?- matrix; 

(2) Fundamental commutation relations for the quantum transfer matrix, involving the 
quantum i?-matrix as "structure constants" , 

which enter the quantum "inverse scattering" method - the first systematic method 
for quantization of completely integrable models. 

In a related development, Drinfeld \W, 13] was the first to realize the deeper al- 
gebraic and geometric nature built-in into the theory of classical and quantum com- 
pletely integrable models. Namely, he showed that the algebraic structures (1) and 
(2), listed above, constitute the basic structural relations of non-commutative and non- 
cocommutative Hopf algebras which ultimately received the name "Quantum Groups" 
and evolved into one of the "hottest" topics in mathematics and theoretical physics. 
Furthermore, it was realized that a quantum group is a deformation of a classical Lie 
group much in the same way quantum mechanics is a deformation of classical Hamil- 



tonian (symplectic) mechanics ||T^. Most importantly, in the "semiclassical" limit the 
basic quantum group algebraic structures (1) and (2), formulated above, transliterate 
into a special distinguished Hamiltonian structure on the classical Lie group G, called 
Lie-Poisson structure, which is compatible with the group multiplication. This is pre- 
cisely the class of Hamiltonian structures given by the quadratic fundamental i?-matrix 
Poisson brackets mentioned above in the context of classical completely integrable mod- 
els. 

The concept of quantum group symmetries in integrable quantum field and sta- 
tistical mechanics' models lead in the recent years to numerous fruitful and exciting 
developments in theoretical physics : from generalization (i.e., q- deformation) of the 
fundamental notions of internal and space-time symmetries in quantum field theory 
and spin-statistics connection to quantum magnetic chains and critical dynamics |[T8[| . 
Furthermore, the concepts of integrability and perturbations around exactly solvable 
theories find their place in (close to) realistic models of elementary particles, e.g. in 
quantum chromodynamics . 

In what follows, few particular topics among those mentioned above are discussed 
in some detail. 

2. Matrix Models of Non-Perturbative Strings and Integrabil- 
ity 

2.1 Conventional Perturbative String Theory 
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The standard geometric formulation of perturbative string theory provides the fol- 
lowing prescription for calculating physical observables (fermionic degrees of freedom 
are discarded for simplicity) : to construct scattering amplitudes one considers func- 
tional integrals over (Euclidean) string world-sheets T^a,g ^ smooth Riemann surfaces 
embedded in D-dimensional (Euclidean) space-time of genus G and area A : 



7 

^string 



Y^g'^ dAe~^^Z^,G (1) 
Za,G = J [Vh] VX exp {-Sstr^ng[X, h]] f{ j d^a,VhV{X, h- h) (2) 



S string ^ I d (7 



^/h{h'^^daX^d,X''G,,{X) + ^{X)R^^\h)) + e'^'daX^d.X^B^^iX) 

(3) 

Here the following notations are used : g denotes the string "coupling" constant, A is 
the "cosmological" constant, the string action Sgtring represents a typical D = 2 confor- 
mally invariant field theory model describing D = 2 gravity (given by the world-sheet 
metric hab{cr) ) coupled to world-sheet "matter" fields X^{a) (describing the embedding 
of T^A,G in )• The functional G^^, ^,B^u represent the space-time dilaton-gravity 
multiplet. String interactions are given in (|l]),(0) by geometrical splitting and join- 
ing of individual string world-sheets (thus creating handles on the total world-sheet), 
whereas asymptotic in-coming and out-going string states are given by vertex operators 
V{X, h; k) {k indicating the momentum of in/out-state). 

Henceforth, for simplicity, we shall suppress the vertex operators Vs in (^, i.e., we 
shall concentrate on the string partition function. 

The most difficult part in calculating (^ is the accurate treatment of the functional 
measure [Dh] over the space of all gauge-inequivalent classes of metrics hab on T,a,g 
w.r.t. reparametrization and Weyl conformal invariance. In case of conformal gauge 
hab = G^hab{T) where (j) is the Weyl conformal factor and hab{T) is a reference metric 
with constant curvature i?*^^-* (Jiab^ and which depends in general on the moduli {r} 
of the corresponding Riemann surface T,a,g , the standard Faddeev-Popov gauge-fixing 
procedure yields ^ P : [Vh] = 5 (hab - e'^habir)) A$n Vcj) (cir) , with the Faddeev- 



Popov determinant A$n giving rise to the well-known conformal anomaly. 

An important result about the entropy of (random) surfaces with fixed area A , first 
obtained by Zamolodchikov [^] in the semi-classical approximation and subsequently 



strengthened in states that for large A : 

Za,g -a^oo constG e^^-^ A'^^'-^^/^^-' (4) 

where Ac denotes a "critical" value of the "cosmological" constant A , x = 2(1 — G) is 
the Euler characteristics of the surfaces, and 70 denotes a critical exponent depending 
on the world-sheet "matter" fields. Relation (H) implies for the string partition function 
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which shows that one can obtain complete nonperturbative result for Zstrmg by taking 
the double scaling limit : 

A ^ A, , g^-^Q such that gl^ ^ [ [K-lf-^^ ^ ^^""^ 



2.2 Lattice Regularization of String Theory; Matrix Model Formulation 

In a series of pioneering papers |^| it was found that statistical mechanical models 



of random matrices ("matrix models" for short) provide an adequate apparatus for 
nonperturbative description of lattice-regularized string theory based on the method of 
random triangulation (and, more generally, random polygonization) of the (Euclidean) 
string world-sheet. A decisive breakthrough occurred further in refs. ||10| which proposed 
ways for correct implementation of the continuum limit as double scaling limit @ 
allowing for exact solutions in string theory. 

Let us note that, whereas matrix model formulation of random surfaces is adequate 
for solving integrable lattice models of planar statistical mechanics, its application to 
genuine string theory is limited so far to the case of Z) < 2 dimensional embedding 
space. Nonetheless, the exact solvability of matrix models provides an important testing 
ground and qualitative hints for the nonpertubative string theory solution in the realistic 
cases (for extensive reviews, see p6[|). 



Since our primary goal here is to elucidate the emergence of integrability structures 
in the matrix models of string theory, we shall consider for illustrative purpose the 
simplest one-matrix model whose partition function is given by : 

Z= d^'M e-^ , V{M) = E ^ Tr (7) 

with M = \\Mij\\ being a N x N hermitian matrix. In ordinary perturbation theory 
defined in terms of a free "propagator" {MijMki) ^q-^ ~ N~^6ikSji and A;-leg vertices with 

weights tkN ij^j , each diagram F gives contribution of the form : 



n I t>^N Uj\ . iv-^(r). iv^(r) = wr [{t}] m^'^ f 



(8) 



A;>3 



On the l.h.s. of (|) Vfc(r) , P(r) , L(r) denote number of k-leg vertices, propagators 
(links) and closed loops (faces) of F , whereas on the r.h.s. x(r) = L{T) — P{T) + 
J2k>3 Va;(F) denotes the Euler characteristics of the two-dimensional polygonized surface 
spanned by F , and Wr [{t}] indicates the product of the vertex weights. Clearly, 
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L(r) = A(r) can be understood as area of F . Thus, the partition function (|^ can be 
written as : 

Z = exp| Y. A^x{r)g-(ln/3/7V)A(r)+ln>Vr[m]l (9) 

liconn. surfaces F J 

i.e., the free energy \nZ of the matrix model (|^) represents the discretized regularized 
partition function of random surfaces (more precisely, "pure" D = 2 gravity with action 
A{r) interacting with "matter" with action InWr ) upon making the following iden- 
tifications (comparing with (|l|)-(0), (|^)) : ~ g ("bare" string coupling constant), 
N/P ~ e-(A-Ac) (^/^ _ "cosmological" constant) and the double scaling limit (^ takes 
the form of a special continuum limit : 



^ — ^ 1 , N — ^ oo , such that g^^^^g 



2-7o' 



1 



fixed (10) 



Explicit solution for the partition function (^ is derived by the method of orthogonal 

polynomials ||25|. Diagonalizing the hermitian matrix M = ?7diag (Ai, . . . ,Xn) and 

1 

integrating over angle variables in U , one gets (rescaling M — > M ) : 

N r ^ 1 

Z= /n^A.A(A)exp -/35:r(A.) A(A) ; V{\) = Y.tuX\, A(A) = J] (A. - A,) 

i=l I i=l J k>0 i<j 

(11) 

Introducing a complete set of orthogonal polynomials Pn(A) = A"+ lower order terms : 

'rfAP„(A)e-'3^(^)P™(A) = hJnm (12) 



and re-expressing the van der Mond determinant as A(A) = (iet||Pj(Aj)|| , the integrals 
in (^) factorize and yield upon using (|T^) : 

N 

Z[{t}]=Uhn{{t}) (13) 

n=0 

where the dependence on the parameters of the random matrix potential is indicated. 
Explicit solutions for /i„ ({t}) can be found through solving the flow equations w.r.t. 
tk which correspond to integrable lattice hierarchies, as briefly discussed in the next 
subsection. 

2.3 Differential Integrable Hierarchies from Matrix Models 

The appearance of integrable hierarchies in the continuum (double scaling) limit ([1 



are extensively discussed in numerous papers (see and refs. therein). It turns 
out, however, that flow equations inherent to integrable hierarchies can be extracted 
directly from discrete matrix models even before taking the continuum limit HE 



which reveals their close connection with topological field theories |29 
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The above result is achieved most easily by employing again the method of or- 
thogonal polynomials (|l^). Namely, on the Hilbert space spanned by {-Pn(A)}„>g one 
introduces two conjugate operators Q, V with matrix elements defined as : 



KQ^n = j rfAP„(A)e-^^WAP^(A) 
KVn.n = J rfAP„(A)e-^^W^P^(A) 



(14) 
(15) 



From (|14D; dUD one easily gets the matrix model "string" equation (the second one 
below) : 



V = P{V'{Q)). 



/3(V^'(Q))(_), Q 



1 



(16) 



where the subscript (— ) denotes taking strictly lower-diagonal part of the corresponding 
matrix. The "string" equation yields recurrsion relations for the matrix elements (jl^) 
of Q . _ 

j) the following fiow equations : 



It is straightforward to deduce from ([T^ and 

dPn{\) 



dtr 



m _ 

dtr 



QnmPmi^) — AP„,(A) 



(17) 
(18) 



Eq. (|T8|) is the integrability condition for eqs. (|T^ and it is compatible with the "string" 
equation (p^. One can easily identify (|T8D with the Lax form of the fiow equations of 
the integrable Toda lattice hierarchy by inserting into (|1^) the explicit form of the Q 
matrix elements : 



Q 



n,n+l 



Qn,n = <Sn-l ; Qn+l,n 



h. 



n+1 



(19) 



the rest being zero as a consequence of the recurrence relations for orthogonal polyno- 
mials. Eqs. ([T8|) are now Hamiltonian, and the lowest one (with r = 1) is generated by 
the well-known Toda lattice Hamiltonian : 



} 



(20) 



Following , it is possible to replace the discrete lattice integrable hierarchy (0) by a 
differential hierarchy at each fixed lattice site n where the continuum variable is x = ti 
. Indeed, from the first (r = 1) fiow eqs. (|1^) and (p!8| ) yielding 



dP, 



n+l 



dti 



TZ P 



dti 



T^n {Sn-1 — Sn) 



dSn 

dti 



(21) 



^ i)n = exp{/rfti5„„i}P„ = hJ-Pn , 

A^n = Qnmi^ra = + Sn-lPn + 7^„_lP„_l) = \d + 7^„ {8 - 5„)~^j (22) 



one obtains (upto gauge transformation P, 
and similarly Q„„ — ^ Qnm = K^Qnmhm ) 
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where d = . Once again using (|2T|), one can rewrite the discrete evolution eqs.(|T^,(|T8|) 
in a differential Lax form for a fixed lattice site n : 

^ = (10.*. . (23) 

where r > 2 , i/j = ipniti) and the subscript + indicates taking the purely differential 
part of the corresponding pseudo-differential Lax operator (cf. last equality in (p2|)) : 



L = d + A{d-B)-' , A = nM , B = SM (24) 



Eqs. (P5D, are immediately recognized as the well-known 2-boson reduction of KP 
integrable hierarchy. 

Using generalization of the method of orthogonal polynomials, it is possible to derive 
flow equations for integrable hierarchies also in the general case of multi-matrix models 
(describing random surfaces interacting with q different types of "matter") : 



/ n^^'M, expi-^ (TrM,M,+i + 5]t,,fcTrMM I 

•' i=l { i=l \ k>0 / J 



(25) 



without passing to the continuum limit ||2^. The appropriate generalization of (|2 
now reads : 

Lg = d + j2Aiid- Bi)-' id - Bi+r)-^ ... (9 - B,)-^ (26) 
1=1 

which is a 2g-boson reduction of KP integrable hierarchy (see subsection 3.3 for more 
details). 

Finally, returning to the string partition function (|13]), one can show that Z [{t}] = 
TN{t) ^ the Toda lattice r-function subject to the so called Virasoro constraints : 

CsZ[{t}] = () , s>-l ; [Cr,Cs] = {r-s)Cr+s (27) 
d " d d d ''''"^ d d 

Otk+s fc=o '^'^s-k Otk-s fc=o '^'^-k Ots+k 

which are equivalent to the constraints on the pertinent integrable hierarchies imposed 
by the "string" equation (|16D and are in fact Ward identities due to the symmetry 5sM = 
EgM'^'^^ , s > — 1 of d^). Similar relations hold for multi-matrix models. Moreover, 
in the continuum (double scaling) limit Z [{t}] can be identified with r-functions of 
reduced KP hierarchies subject to the so-called l^-constraints (generalizations of (|27|) 
which span VT-algebras; cf. (|iD|)-(p(]|) below). 



3. Integrable Systems in Classical Physics : Geometric Formu- 
lation 
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Now we pass to a brief review of some principal structures and properties of integrable 
systems in classical and quantum mechanics (or field theory) . We first recall the notion 
of integrability. 

Complete integrability: Consider a Hamiltonian system with n degrees of freedom pos- 
sessing standard Hamiltonian structure with a Hamiltonian H{p, q) and Poisson bracket 
{■,■}. A Hamiltonian system is called completely (or Liouville) integrable if it has n in- 
dependent conserved quantities (integrals of motion) /^(p, g) , = 1, . . . , n, which are 
in involution: {Ii,Ij} = 0. For such a system we can find the action-angle canonical 
variables and write explicitly the general solution to the equations of motion. 

Lax formulation: For infinite-dimensional (field theory) integrable Hamiltonian systems 
, there exists the convenient Lax (or "zero-curvature") formulation In the Lax 
formulation, the phase space of the Hamiltonian system is parametrized by elements 
L taking values in some Lie algebra Q and the dynamical equations of motion can be 
written in terms of a Lax pair L, P, the latter similarly taking values in Q, as the 
Lax-type equation 

f^lL.P] (28) 

The Lax formulation leads straightforwardly to the construction of the involutive inte- 
grals of motion. Namely, for any Ad-invariant function I on Q, I (L) is a constant of 
motion. In fact, it can be shown that any completely integrable Hamiltonian system 
admits a Lax representation (at least locally) . 



3.1 Adler-Kostant-Symes/Reyman-Semenov-Tyan-Shansky Scheme 

A very wide class of integrable models can be constructed through the application of 



the AKS-RS approach ||3T|, [T5| having roots in the group coadjoint orbit method [p^] . 

Let G denote a Lie group and Q be its Lie algebra. G acts on Q by the adjoint action: 
Ad{g) X = gXg~^, with g & G and X & Q. Let Q* be the dual space of Q relative to a 
non-degenerate bilinear form (■!•) on Q* x Q. The corresponding coadjoint action of G on 
g* is obtained from the duality of (•!•): {Ad*{g)U\X) = {U\Ad{g-^)X). We will denote 
the infinitesimal versions of adjoint and coadjoint transformations (for g = expY) by 
ad{Y)X = [Y,X] and {ad*{Y)U\ X) = - {U\ [Y , X]) , respectively. In particular, 
when Q is endowed with an ad-invariant bilinear (Killing) form (-, ■) allowing to identify 
g* with g , we have ad*{Y)U = [F, f/] . 

There exists a natural Poisson structure on the space (^*, R) of smooth, real- 
valued functions on g* (sometimes called Kirillov-Kostant (KK) bracket), given by : 



{F,H}{U) = -{U 



VF{U) , VH{U) ) (29) 



where F,H e (^*,]R) , the gradient VF : g* — > g is defined by the standard 
formula -^F {U + tV) \t=o = {V \ VF([/)) and [■ , ■] is the standard Lie commutator on 
g. On each orbit of G in g* the Poisson bracket (^) gives rise to a non-degenerate 
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symplectic structure. Moreover, for any Hamiltonian function H on such an orbit we 
have a Hamiltonian equation of motion dU/dt = ad*{VH{U))U ( = ['S/H{U) , U ] when 
Q has Kilhng form) . 

Ref.[p!5[ introduced the -R-operator (generahzed i?-matrix) as a hnear map from a 
Lie algebra Q to itself such that the bracket: 



[X,Y]n^\[RX,Y] + ]^[X,RY] 



(30) 



defines a second Lie-bracket structure on ^ , or equivalently, defines a second Lie algebra 
Qn isomorphic to ^ as a vector space. The Jacobi identity for the /^-commutator ( pO| ) 
implies that the modified Yang-Baxter equation (YBE) for the i?-matrix must hold (for 
arbitrary Xi_2,3 ^ Q ) '■ 



E 

cyclic(l,2,3) 



Xi , RX2 , -RX3 



R 



RX2 , X3 



+ 



X2 , -RX3 







(31) 



A sufficient condition (occurring in all interesting cases) for the fulfilment of (|3TD is : 



RY 



r] + [x, RY 



—a 



X, Y 



(32) 



[a being arbitrary constant), which usually is written in terms of the "ordinary" r- 
matrix r & Q ®Q canonically isomorphic (upto a factor) to R E Q ®Q* via the Killing 
form : 

(33) 



r(i2) 


(13)1 




r(i3) 


(23)1 




r(i2) 


(23)] 


r , 


r 


+ 


r , 


r 


+ 


r , 


r 



(12) 

r E 



a(i— invariant 
r,,T ®leU{g)®U{g)®U {g) etc. 



where lA {g) is the universal enveloping algebra. Classification of solutions of (|33|) 
(depending on a spectral parameter) for simple g is given in refs. ||33|| . 

With the help of (^) one can furthermore introduce on ^ ^* a new KK-type 
Poisson bracket {■,-}r called i?-bracket by substituting the usual Lie commutator [-, ■] 
for the i?-Lie commutator [■,-]r ( pOD in (p9|): 



{F,H}n{U) = ~{U 



VF{U) , VH{U) 



R 



(34) 



A function H on ^* is called 74d*-invariant (Casimir) if H[Ad*{g)U] = H[U] or, 
infinitesimally, ad* {V H {U)){U) = for each U E g* . Then one can show [|T5[ that: 



(1) the a(i*-invariant functions are in involution with respect to both brackets (29) and 



(2) the Hamiltonian equation on g* ~ g*^ takes the following (generalized Lax) form: 

]^R {VH{U)) , f/] for ^ with Kilhng form^ 
^ (35) 



dU/dt= -ad*[R{yH{U)))U 



= TyT* ® and \RX ^ T'r.,j {T^,X) , where {T*} denotes a basis in G 
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Eq. (|35| ) can be obtained from variational principle with the following geometric action: 



dU 



W[U] = - 
ad*ji {yR{U))U 



U 



dyR = 



dtH[U] 



(36) 
(37) 



where the integrals in (|36| ) are along arbitrary smooth curve on the phase space Q* ; 
H[U] is any Casimir on Q* ; yR{U) is Maurer-Cartan one-form on Qr and function of 
U E Q* determined by the first eg. (|37|) , with the -R-coadjoint action : 



(3J 



ad*^{X)U = -{ad\RX)U + R*{ad*{X)U) 



2 

RX, U 



X, U 



for Q with Killing form 



Hence the above AKS-RS technique leads to a direct construction of completely 
integrable systems where the set of independent Casimir functions on Q* is precisely the 
complete set of integrals of motion in involution. The basic realization of this technique 
arises when the Lie algebra Q decomposes as a vector space into two subalgebras 
and i.e. Q = Q) G~. Let P± be the corresponding projections on Q±. Then 
R = P+ — satisfies the modified YBE (|3TD and provides a specific realization for the 

), (|35|) and (|38|) take the following simple form: 



above scheme. In this case eqs.(pOD, 



X,Y 



X^, 



{ad*j,{X)U 



± 



x+,u^ 



dU 



= T 
'5H 

Ju 



± 



(39) 
(40) 



where x± = p±x e g± , u± = p;u G (0^)* , [X, u]^ = p;[[x , U]) e {Q: 



3.2 Algebra of Pseudo-Differential Operators and Integrable Hierarchies of 
Kadomtsev-Petviashvili type 

Here we will illustrate the AKS-RS construction on ^ = '^T>0 - the algebra of pseudo- 
differential operators on a circle. Recall that an arbitrary pseudo-differential operator 
X(x, Dj;) = J2k>-oo ^k{x)D^ is conveniently represented by its symbol |3^ - a Laurent 
series in the variable ^ : X(^,x) = J2k>-oo ^k{x)^'' , and the operator multiphcation 
corresponds to the following symbol multiplication: 



x(e,x)oy(e,x) = J2 



N>0 



1 d^Xd^Y 



(41) 



which determines a Lie algebra structure given by : [X, Y] = (X oY — Yo X). 
On ^)T>0 one introduces an invariant, non-degenerate bilinear form: 

L\X) = Tta (LX) = J dx Res^ ( L(^, x) o X(^, x) 



(42) 
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which allows identification of the dual space ^T>0* with ^!T>0 . 

There exist three natural decompositions of ^ = ^T>0 into a linear sum of two 
subalgebras Q = Q^j^® Qi_ labelled by the index £ taking values £ = 0, 1, 2: 

{oo I oo 

i=e ) y i=-l+l 

(43) 

Correspondingly the dual spaces to subalgebras are given by: 

oo I I oo 

L_ = L<_, = J2 D-'o u.,{x) \ ; gi* = \L+ = L>_, = ^ ° ^i(^) 

i=£+l J [ i=-e 

(44) 

Note that in ( ^^ the differential operators are put to the left. Henceforth, we shall 
skip the sign o in symbol products for brevity. 

Defining Ri = P+ — P_ for each of the three cases, eqs. (BUf) take the form : 

[X, Y]r^ = [X>e, Y>i] — [X^i, Y^i] , ad*j^^{X)L = [X>i, L^^i]^_^ — [X<^, -^^>-^]>_^ 

"(45) 

Now, choosing an infinite set of independent Casimir functions : 

H^+i = I dx Res L^+i , m = 0, 1, 2, . . . (46) 

m + 1 J 

the three decompositions ( ^i)) of ^VO labelled by £ = 0, 1, 2 yield, according to the 
AKS-RS scheme, three different integrable hierarchies - the standard KP hierarchy 
{I = 0) and the first and second modified KP hierarchies {i =1,2) : 



dL 

dtm. 



+ 



(L'")>,,L =0 (47) 



One can show that all three KP hierarchies are related through symplectic gauge trans- 
formations [^. Here we shall concentrate on various Poisson reductions^oi the stan- 
dard KP hierarchy which appear, in particular, in the context of matrix models of 
strings. 

The phase space of the standard KP integrable system is : 

Mkp = |l = D + f; n,(x)D-^j (48) 

''Let us recall the general notions |Q . Let {M. ,P) be a smooth Poisson manifold with Poisson 
structure P : T*{M) — > T{Ai) and let 5 be a smooth submanifold of A4 with embedding 
^ : S — > M . Now, a Poisson structure P' : T*{S) — > T{S) on S C M is called Poisson 
reduction of P if for arbitrary functions on M the following property is satisfied : fi* ({/i, /2}p) = 
/i*/2}p' • In other words, restriction of the Poisson brackets w.r.t. P of arbitrary functions on 
M to the submanifold S is equivalent to computing the Poisson brackets w.r.t. P' of the restrictions 
on 5* of these same functions. 
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(it by itself is a trivial Poisson reduction from Ai = '^T>0* ^ ^VO) . The KK Poisson 
brackets (first eq. (^OD ) read on (^) : 

{Un{x) , Um{y)] = Qn-l,m-l{u{x))S{x - y) (49) 
nnmHx)) = J2(-'^)[l^j^n+m-k+lix)D^ ~^[k j ^^^n+m-k+l{x) (50) 

and are recognized as the (centerless) Wi+oo algebra (see 0), which is isomorphic to 
VOV{S^) C ^PO (algebra of purely differential operators) 

Now, returning to the flow eqs. in (multi)matrix models ( pSj) with Lax operator 



6]), one can show ||3^ that the space of 2g-boson Lax operators : 

M,= \L, = D + j2Ai{D- Bi)-' {D - Bi+,)-' . . . {D - B,)-'] (51) 



1=1 



is a legitimate Poisson reduction of the full KP hierarchy given by (|^), i.e., the flow 
eqs. (p3D indeed correspond to Hamiltonian integrable systems. Furthermore, what is 
even more important, one can express all 2q fields parametrizing (|5TD in terms of Dar- 
boux canonical pairs of coordinates (^ar{x), br{x) 



1 

r=l 



Bi = bi + bi+i + ---bg , Ag = ag (52) 

g— 1 »13— 1 n2 — 1 

Aq.ria, &) = E ■ ■ ■ £ E + K + ■■■+ bnr-r+l) " " " (<9 + &ni ) (53) 

nr=r n2=2ni=l 

{ar{x), bs{y)]p, = -6rsdJ{x - y) (54) 
Expressing Lq in (|5ll) as power series in : 

oo 

Lq = D + Y.Uk[{a,b)]{x)D~'' (55) 

k=l 

min(g— 1,A:— 1) q 

Uu[{a, b)]{x) = aqP^'\ [bg) + b)Pk-i-rK \ + Vi, • • • , E ^0 (56) 

r=l l=q—r 

where Aq_j.{a^b) are the same as in (|53|), and -P^^-* denote the (multiple) Faa di 
Bruno polynomials : 

Pi^)(i?iv, i?iv-i, . . . , i?i) = E (-^ + ^i)"' ■ ■ ■ (-5 + BmT^ ■ 1 (57) 

miH ^m,ff=n 

one obtains a series of explicit (Poisson bracket) realizations of Wi+oo algebra in terms 
of finite number of 2q bosons (a^, K)l=i for any g = 1, 2, 3, . . . (cf. (^)-(|50|)) : 

{Uk[{a, b)]{x) , Ui[{a, b)]{y)} = Qk-i,i-i{U[{a, b)Mx - y) (58) 
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Similar "Darbouxization" of Lq (^) holds w.r.t. the second KP Hamiltonian structure 
(see next subsection). 

3.3 Lie-Poisson Groups and Lie Bi- Algebras : Hamiltonian approach 

Completely integrable systems possess another notorious feature - that of multi-Hamil- 
tonian structure, i.e., they always possess at least two independent compatible Poisson 
structures In the case of KP integrable hierarchies (and their reductions) it turns 
out that the "second" Poisson structure can be understood in terms of the Lie-Poisson 
structure on the Lie group G corresponding to the Lie algebra Q entering the AKS-RS 
construction, whereas the "first" Poisson structure is just the KK structure ( PH) on 

The notions of a Lie-Poisson group and its "infinitesimal" version - Lie Bi-algebra 
were first introduced by Drinfeld |jT6|. A Lie group G is called Lie-Poisson if on the alge- 
bra of smooth functions on it Fun{G) there exists Poisson structure which is compatible 
with the group multiphcation, i.e., A (^{ Fi , F2 }) = { A (Fi) , A (F2) } where A de- 
notes the coproduct in Fun{G) ( AF{g, h) = F{gh) ^. More precisely, the Lie-Poisson 
structure is given by : 



{ F,{g) , F^ig) } = (vLF^ig) ® VLF2{g) 



r{g) 



^LF{g) 



t=o 



^RFi{g)®VRF2{g) 
^RF^g) 



i=0 



Fige 



(59) 
(60) 



where Vl,r denote left/right Lie-derivatives, and where r{g) is a cocycle on G with 
values in Q : 

r{gh) =r{g) + Ad{g) ^ Ad{g)r{h) (61) 
In the case when r{g) is a cohoundary , i.e. 



r{g) = Ad{g) O Ad{g) tq - tq 



(62) 



with To e Q ^ Q being a constant element, one easily finds that the Jacobi identity for 
(|59D reduces precisely to the YBE (^) for tq as a classical r-matrix. 

Let us note that for matrix groups eq. (|5UD can be written in a simpler form : 



{gfg} = r{g)g®g 



{g^g} 



1 g ® g] in case of (62) 



Eq.(^) implies the following exterior derivative equation 



dr{g) 



[Y{g)®l + l®Y{g),r{g)]-(t){Y{g)) 
d 



0(A) 



dt 



r 

t=o 



(e^*) for eg 



(63) 

(64) 
(65) 



^Two Poisson structures {•, - ji and {•, •}2 are called compatible if their linear combination {•, -Jx 
{■, -ji + A{-, •}2 is also a Poisson structure. 
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where Y{g) = dg g ^ denotes the Maurer-Cartan one-form on Q , and 0(-) defined in 
(|65D is a ^ ® ^-valued co cycle on Q : 

4>{[x ,Y]) = [x ®i + i®x , 4>{Y)] + , r ® 1 + 1 ® r] (66) 

The solution of (Q) is a coboundary r{g) ( |62D if and only if 0(-) is coboundary : 

0(X) = [ro , X ® 11 + 11 ® X] (67) 

The cocycle 0(-) allows to introduce a Lie-commutator [■, -j^ on the dual space Q* and, 
correspondingly, coadjoint action adl{-) of ^* on Q as follows : 



t/, V 



X 



V 



ad:{U)x) 



U0V 



0(X) 



yu,v e g* , \fx e g (68) 



whereas the mixed commutator between elements of g and g* is defined as : 

[X , U] = ad*{X)U - adl{U)X (69) 

According to an important theorem by Drinfeld |T^, the Lie algebra g of each Lie- 
Poisson group G possesses a Lie-bialgebraic structure and vice versa. Namely, by means 
of ( p8| ) and ( |6PD the direct sum (as vector space) V = g(Bg* becomes itself a Lie algebra 
called "the double" , such that g and g* are isotropic subalgebras of V w.r.t. the Killing 
form on V : 



{{X,,U,) , (X2,f/2)) = (f/ll X2) + (f/21 Xi) V(Xi,2,?7i,2) G V 



(70) 



The triple (^P, g, g*^ is also called Manin triple. Furthermore, the double group D ^ 
G X G* , corresponding to P , is direct product (as a manifold) of its subgroups G and 
G* corresponding to g and g*, respectively. Let us stress the symmetry of the above 
construction under the exchanges g < — > g* , G < — > G* . 

Now it is easy to write down the explicit solution for the Lie-Poisson group cocycle 
rig) (ill) : 



{U^V\r{g)) = Ug-'Vg) 



g-'ug 



Vf/, V eg* 



(71) 



where the subscripts (±) indicate projections in the double algebra V along g , g* , 
respectively. 

As in the AKS-RS Lie-algebraic scheme, one can construct integrable Hamiltonian 
systems on Lie-Poisson groups. Namely, from (31) or (|62D one can easily verify that 
all y4(i(-)-invariant functions on G , H[hgh^^] = H[g] , are in involution w.r.t. the 
Lie-Poisson structure (^) : {Hk[g], Hi[g]} = . Similarly, the analogues of the 
Hamiltonian eqs. of motion ( |53| ) and the associated geometric action (^) take now the 
form : 

dg 



dt, 



-9 



s.iYig)) 



y{9) 



dtHk[g] 



(72) 
(73) 



String Theory and Integrable Systems 



16 



where the action of the operator fg : Q* — > Q is defined by (f/| fg {V)) = {U ®V\ r{g)) 
cf. (fflD , Sg = f~^, and denotes inverse operator of the exterior derivative acting on 
the corresponding closed 2-form ^ (recall Y{g) = dgg^^ ). 

The principal example is provided by the extended Volterra algebra of purely pseudo- 
differential operators (c is arbitrary constant) : 



g = {^VO)_ = \j2Mx)D-'' + c\nD 

U>i 



(74) 

whose Lie-double is the extended algebra of all pseudo-differential operators : 



l>0 



^ (75) 

Here E indicates the central element of Wi+oo , as well as of the whole ^VO , which is 
dual to In D , cf . [0, ^ . The corresponding extended Volterra group (exponentiation 



of (17^., 

G = (vroo) _ = K = L= il + Y^ Uk{x)D-^ ) o d4 (76) 



k>l 



can be viewed as a set of spaces (for each c fixed) of Lax operators of generalized 
KP hierarchies w.r.t. the second KP Hamiltonian structure |^ (eg. (|77D below). The 
latter is given precisely by the Lie-Poisson structure (|63|) with the cocycle r{g) ( ffl]) for 
G = (© (see 

Let us go back to the example of 2g-boson KP Lax operators appearing in the multi- 
matrix string models (pl|). The second Hamiltonian structure for general Lax operators 
has the formf] : 



{ {L\ X) , (L| Y) } = Tr^ [{LX)^ LY - (XL)^ Y L) + 

y"rfxRes([L, X])9"^Res([L, F]) (77) 



Here Tr^ denotes the Adler trace (|i2| ) and the subscript + indicates taking the purely 
differential part. For the coefficient fields Uk{x) of L (|48|) the second KP Poisson algebra 
( iTTf ) yields the nonlinear (i.e., non-Lie) Wqo algebra The latter appears as a unique 
(modulo certain homogenity assumptions) nonlinear deformation of Wi+oo algebra. 

In analogy with eqs.(^2|)-(^) one can express ^l|] the coefficient fields {Ai,Bi)l^^ 
of Lq ( pID in terms of Darboux canonical pairs of fields (c^, GrYr^i w.r.t. the second KP 

^Closeness follows from the exterior derivative equation dsg = ad*(F((7))sg — Sgad(y(g)) — 
Sg4'{Y{g)')sg , satisfied by the inverse cocycle operator Sg as a consequence of eq.(|6^). 

^The second term on the r.h.s. of (|7^) is a Dirac bracket term due to the second class constraint 
1*0 = in L p 
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Hamiltonian structure ([77| ) : 

5fc = Cfc + ^ Q , l<k<q ; Ag = J2 (9 + Cr) Cr (78) 

l=k r=l 
k / rik+q-k 

^k= J2 \9 + en,- en,+g-k + Q, 
"fc=l \ lk='n.k 
rik / rik-i+g-l-k \ 

X £ id + en,_, - en,_,+g-i-k + E c,,_Jx--- 

nfc_i=l \ 'fc-i=rafe-i / 

ns ?12 

X H (9 + e„2 - e„2+i + c„2 + c„2+i) 51 ("^ + ^"1 ) 1 < ^ < 9 - 1 (79) 

n2=l ni=l 

{cfc(x) , e;(?/)} = -5m9^5(x - ?/) , fc,/ = 1,2, . . . ,g (80) 

These equations are equivalent to the following "dressing" form for the 2g-boson KP 
Lax operator : 

L, = D + j2Ai{D- Bi)-' {D - Bi+,)-' ...{D- B,)-' = Ug . . Mi DV{\ . . V^' (81) 
1=1 

Uk = {D-ek)eS"' , Vfe^e/^M^-efe) , fc = l,...,g (82) 



Eqs.([78D-(^) or, equivalently, eqs. (|8TD-(^2D can be viewed as generalized Miura trans- 
formation for the 2g-boson KP hierarchy 

The Miura-transformed form of Lq (|8lD reads explicitly : 



Lq = D + Y^Uk[{c,e)]{x)D-^ (83) 

k=l 

f/4(c, e)]{x) = (e, + c,) ^2(9 + q) d + 



min(5— l,fc— 1) q 

Aq_r{c, e)Pjf^^^X{eq + Cg, eg_l + Cg_l + Cg, . . . , Cq-r + I] q) (84) 
r=l l=q—r 

where Ag_r(c, e) are the same as in ([79D, and -P^^-* denote the (multiple) Faa di 
Bruno polynomials (0). 

Now, in complete analogy with eqs.(^)-(|58|), which yield a series of realizations 
of the linear Wi+oo algebra in terms of 2q bosons, we obtain, upon substitution of 
(|83D-(^) into ([77|), a series of explicit (Poisson bracket) realizations of the nonlinear 
Woo algebra in terms of 2q bosonic fields for any g = 1, 2, . . . . This algebra plays an 
important role as a "hidden" symmetry algebra in string-theory-inspired models with 



black hole solutions 42 



®For discussion of the generalized Miura transformation and the associated Kuperschmidt- Wilson 
theorem, we refer to [Eol. 
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Concluding this section, let us note that in the general case Lie-Poisson groups pro- 
vide natural geometric description of the dressing symmetries in completely integrable 
models H^, WM . Another outstanding role played by Lie-Poisson structures (Bl) is their 



appearance in the context of the classical inverse scattering method as fundamental 
Sklyanin brackets for the monodromy matrix g ~ T(A) of the auxiliary linear spectral 
problem. 

4. Quantum Integrable Models 

4.1 Quantization of Lie-Poisson Groups : Quantum Groups 

Historically, quantization of completely integrable models, whose Hamiltonian structure 
is based on group Lie-Poisson structures, laid for the first time to explicit construction, 
in the context of the quantum version of the inverse scattering method ||Tl| , of quantum 



groups which were subsequently identified with quasi-triangular Hopf algebras |T3|. 

Among the various ways to introduce quantum groups there exists an approach |12|, 
whose conceptual point of view underscores both the quantum mechanical as well as 
the Hopf algebraic aspects in quantization of Lie-Poisson groups. Namely, on one hand 
Fun{G) can be viewed as Abelian associative algebra of "observables" of a classical 
Hamiltonian system {A4,V) with a phase space M. = G and Poisson structure V = Vlp 
given by (^) : 

Vlp F2) = { Fi , F2 } = (VlFi ® VlFs - V rF^ ® V^Fal tq) (85) 



where the classical r-matrix tq satisfies the classical YBE (|33|) (from now on we shall con- 
sider only coboundary Lie-Poisson groups (|62D). On the other hand, one can easily check 
that Fun{G) is endowed with a structure of a commutative, but non-cocommutative, 
Hopf algebra Ao{m, A, S,e) = Fun{G) with a product m{Fi,F2){g) = Fi{g)F2{g) , 
coproduct A{F){gi, g2) = F{gig2) , antipode SF{g) = F{g~^) and counit e{F) = F{e) , 



and this Hopf structure is compatible with the Poisson structure (|85D , i.e., AoV = VoA 



Thus, quantization of a Lie-Poisson group G may be viewed as a generalization of 
Weyl quantization Fun{G) — > Funh{G) of a classical Hamiltonian system defined by 
(Ai = G,Vlp) , -i.e., non- commutative deformation of the product m(-, ■) — >• m/i(-, •) 
with a deformation parameter h, which satisfies the additional condition that the 
deformed algebra FurihiG) = Ah{rnh, A, S,e) is again (non-commutative and non- 
cocommutative) Hopf algebra which is a deformation of Ao{m, A, S*, e) = Fun{G) . 

Let us recall [|17|, that in ordinary Hamiltonian mechanics on a Poisson manifold 
{Ai,V) with local coordinates (x*) and constant Poisson tensor { , x-' } = P*-' , Weyl 
quantization is given by the associative and non-commutative Moyal product : 

nif, (Fi, F2) = F,^hF2 = mo e^^ (Fi, F2) = F, ■ F2 + ^ { F, , F2} + 0{h'') (86) 

V = P'^TT- ® -TT- ■■ Fun{M) ® Fun{M) — ^ Fun{M) ® Fun{M) (87) 
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where V is the Poisson bi-vector field. Recall also, that the form of the first order 
term in the /i-expansion (last eq.(p6D) is dictated by the semidassical correspondence 
principle. 

In the case of Lie- Poisson groups (G, Vlp) , the deformed product mh{-, ■) preserving 
the Hopf algebra structure and satisfying the semidassical condition, can be constructed 
as follows [0. Let us choose a basis {X^} in U (Q) - the universal enveloping algebra 
of the Lie algebra Q of G, and let hl^r denote the representations of U (Q) in terms of 
left/right Lie derivatives : 7?^^/? (X*) = V^^ (see eq.(pO|)). Then : 

m,,(-,-) = moA , A = (vtl ® vtl) (A) o (vr^ ® ttk) (A^^) (88) 

dimQ dimQ „ l 

A(X,r)= Y: C{.},mih) n {^T n {r^r = t+-r,,X^Y^ + 0{h') (89) 

{a},{/3} i=l j=l 



with the following notations. The coefficients in A(-, ■) : U {Q) ®IA {Q) — > U {Q) ® 
U (^) [[/;,]] are power series in h ; {X*} and {l^-*} are generator basises in the first and 
second copy of lA (^), respectively; ||rjj|| = r^ is just the classical r-matrix satisfying 
(|55D. Moreover, the associativity condition for m/i(-, ■) ( PB[ ) implies the following basic 
quadratic equation on U. [ff) ®U {Q) ®IA {Q) [[/;,]] ( X, F, Z below correspond to the first, 
second and third factor U {Q) in the tensor product) : 

A(X + y,Z)A(X,F) = A(X,F + Z)A(F,Z) ; A(X, 0) = A(0, F) = 1 (90) 

Defining i?(X, Y) = A-\Y, X)A(X, Y)eU{g)®U (<3) [[h]] , one obtains from (H) : 



R{X, Y)R{X, Z)R(Y, Z) = R{Y, Z)R{X, Z)R{X, Y) (91) 
R{X, Y)R{Y, X) = l ; R{X, Y) = l + h r^^X'Y^ + Oih^) 

R is called universal quantum R-matrix associated with the classical r-matrix \\rij\\ = tq 
. If p : Q — > Endiy) is some representation of ^ in a (finite-dimensional) vector 
space V , then the matrix R= {p® p) (^R^ G End{V ® V) satisfies the famous quantum 
Yang-Baxter equation (QYBE) (plus the "unitarity" condition, P G EndiV ® V) being 
the permutation operator) : 

(12)(13)(23) (23)(13)(12) 

RRR = RRR ; RPRP = 1 ; R = 1 + hro + 0{h^) (92) 

The indices (12), (13), (23) indicate the various embeddings of i? G EndiV ® V) in 
End{V ®V®V). 

In particular, from the deformed product (p8D mh (i^i, F2) = m o A [Fi, F2) where 
the functions ^1(5') = Qab , -^2(5') = dcd are just matrix elements of the group element 
g E G (for matrix groups), one obtains using matrix tensor notations : 



R{g ® 11)(11 ® c/) = (1 ® g){g ® 1)R 



(93) 
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The semiclassical limit of ( p3D is precisely given by the Lie-Poisson bracket (|63D. Eq.( 
is nothing but the famous fundamental commutation relations for the matrix elements 
of the quantum monodromy matrix (with spectral parameter dependence suppressed) 
in the quantum inverse scattering method ||TT|. 



For various treatments and numerous applications of QYBE, see ref . |^5| . For parallel 



developments in the abstract Hopf algebraic context, we refer to [^ |. 

4.2 "Soliton" Scattering in Completely Integrable Models 

Let us now consider integrable D = 2 relativistic field theories whose Hamiltonian dy- 
namics is given by actions S[(j)] = J d'^x C{(f), dcj)) - local functionals of the fundamental 
fields (collectively denoted by (t){x^,x~) ) and their derivatives. As usual, one uses the 
light-cone form of the space-time coordinates : = | (x^ ± x°) . Complete integrabil- 
ity in this context implies the existence of an infinite number of independent integrals 
of motion in involution Q^^'^ , whose densities are local (as functionals of and its 
derivatives) conserved currents : 

= j (r^'+^^dx^ + e(^-i)dx+) , s = 1, 2, . . . ; a+T("+i) = d.e^'-^^ (94) 

(here s indicates the D = 2 Lorentz weight). Thus, quantization of completely inte- 
grable field theories means fulfilment of the quantum renormalized Ward identities for 
the renormalized quantum conserved currents (t^^~^^\ Q^'^'^^^ : 

d+ (T^'+^\x)(j){xi) . . . (t){xrS)-d- (Q^'~^\x)(t){xi) . . . 0(x„)) = (5-function terms (95) 

(as usual, (. . .) here denote time-ordered correlation functions). The infinite set of Ward 
identities lead to severe restrictions on the particle ("soliton") scattering processes - 
conservation of all (odd) powers of momenta of incoming and outgoing particles : 

Nin Nout 

1=1 1=1 



which, in turn, implies |^ : (a) no multi-particle production, i.e., Nin = Nout ] and (b) 
factorization of multi-particle scattering amplitudes. The latter property is of tremen- 
dous importance, as it leads to the remarkable Zamolodchikov's factorization eqs. for 
the 3-particle amplitudes 0, meaning that any 3-particle scattering process is accom- 
plished as a sequence of 2-particle scatterings only and, moreover, the amplitude does 
not depend on the order in which these sequential 2-particle scatterings occur: 

^hi2 {^12) '^fciia ^^isj Ok2k3 1^23 j - -^ijig [U23) ^i^ks l^lSj ^kik2 1^12J {^'> 

with 9ab = Oa — Of, , a, 6 = 1, 2, 3 . In (0) the following notations are used : 5*^^' (6*12) 
denotes 2-particle scattering amplitude of incoming particles of "type" labelled by the 
indices i and j and with (on-mass-shell) momenta pi^2 = ""^1,2 (cosh6'i^2, sinh6'i^2), re- 
spectively ( ^i_2 are the relativistic "rapidities"). 
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Now, denoting by V the vector space of internal particle symmetry (particle "types" ), 
one can regard the matrix of the 2-particle amplitude as : 

S (612) = ll^f/ (^12) II e Mat{V) ® Mat{V) (for fixed Ou) (98) 

(12) 

S (^12) = S (^12) ® 1 G Mat(V) ® Mat(V) ® Mat(V) 

(13) (23) 

and, accordingly, for S (6*13) and S* (^'23) • Then it is straightforward to identify (|97|) 
with the QYBE ( p2D in the quantum group framework. 

A closely related natural appearance of quantum group structure in "soliton" scat- 
tering is provided by the notion of asymptotic states' symmetry |^ : \{0,i)) 



TijiO)\{ej)) ; \{e,i)) e V , ||Ty(^)|| e Mat{V) , as a result of the integrability. 
Namely, the 2-particle S'-matrix : 

I (^1,^1) , (^2,^2))^" = sii^^ (^12) I (^2,j2) , ie„j,)r' (99) 

(12) (1) (2) (2) (1) 

can be viewed as a mapping (for fixed "rapidities") S '■ V V — ® V and, 
therefore, its invariance under the asymptotic states's symmetry : 

(12) (1) (2) (2) (1) (12) 

S (^12) T (^1) T (^2) =T {62) T (e,) S (eu) (100) 



is straightforwardly identified with the structural relations ( p3D for quantum groups. 

Let us point out that quantum group relations of exactly the same form as ( P?]) and 
( |100|) do appear in exactly solvable lattice models of planar statistical mechanics, but in 



this case - with purely imaginary "rapidities" { 9 = ia , a being angles characterizing 
the rectangular lattices), S^j (a) being the matrix of Boltzmann weights at each lattice 
vertex, and Tjj(a) denoting the row transfer matrix 



4.3 Quantum Field Theory Approach to Integrable Models with Dynamically 
Broken Conformal Invariance 

Finally, let us briefiy discuss construction of higher local quantum conserved currents 
fulfilling the Ward identities (^), which is the heart of the quantum field theory ap- 
proach to quantization of completely integrable models. Recently, Zamolodchikov [p^ 
proposed powerful general formalism based on treating integrable models as mass per- 
turbations of conformal field theories : S[(f)] = Sconf[(p] + J2iJi^i I d'^x Bi{(f),d(j)) , where 
the "coupling" constants have positive mass dimensions and -Bi(0, (90) are composite 
fields with conformal dimensions less than 2. Since in general it is not possible to find 
explicit expressions for S'con/[0] and Bi{(f), d(f)) as local functionals of local fundamental 
fields {4>}, Zamolodchikov's approach is purely algebraic using results from represen- 
tation theory of Virasoro algebra, in particular, information about the spectrum of 
conformal field dimensions. 
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There exist, however, interesting classes oi D = 2 integrable field theories, i.e., 0{N) 
nonlinear sigma-models and their supersymmetric generalizations with Lagrangians : 



NLa 



n 



N/g 



n 



n 



N/g 



(n\ . . . n^)(101) 
, = (1102) 



which are conformally invariant on the classical level but, upon quantization, they un- 
dergo dynamical dimensional transmutation, manifested through dynamical mass gen- 
eration, leading to anomalous conformal symmetry breakdown. This clearly precludes 
the use of conformal perturbation approach to ( |101| ), ( |102| ). Fortunately, there exist 
long ago alternative nonperturbative f\ treatment of quantum field theory models with 
0{N) or SU{N) internal symmetry - the expansion |51 . 

Let us briefiy illustrate the construction |Q of higher local quantum conserved 
currents for ( [L01| ) [^] within the expansion framework (the same techniques applies 
to other 1/A^-expandable integrable models as well). The 1/N expansion is obtained 
from the generating functional of time-ordered correlation functions : 

Z[J] = Jvnl[6{n^- N/g) exp |z j (fx ]^{dnf + (j, n) 



Va exp |-yS'i[o-] 
SM = Trln( 



(fx a 



' 2 

-92 + 0- 



(103) 
(104) 



by expanding the effective a-field action (|104|) around its stationary point cTc = m = 
,2„-in/g ((;iynamically generated mass of the "Goldstone" field n , fi being the renormal 



/i e 

ization scale), i.e., a{x) = rn^ + --^a{x) . As a result, one arrives at the 1/A^ diagram 
technique with (free) propagators in momentum space : 



(0) 



/ 2 2 
[m + p 



5' 



ab 



(^^)(0) 

1 



(105) 



(27r)2 (m2 + F) (m2 + (p-A;)2) 



and ordinary tri-linear arm-vertices. 

It turns out |^4[, that the 1/A^ expansion can be renormalized by adapting the well- 
known BPHZ^ renormalization technique. Another remarkable property of the 1/A^ 

^The term "nonperturbative" refers to expansions different from ( or, e .g., pa rtial resummations of) 
the ordinary perturbation theory w.r.t. the coupHng constant g in ( 101 ) and ( 102 ) which is plagued 



by infrared divergences in Z? = 2. 

Bogoliubo v-Pa r asiuk- Hepp - Zimmermann As shown in [Q, the (supersymmetric) nonlinear 
sigma-models (101) and (|l02| ) are renormalizablc within the 1/N expansion also in D = 3 space- 
time dimensions in spite of their naive nonrcnormalizability w.r.t. the ordinary coupling constant 
perturbation theory (note e.g. the presence of the four-fermion term in (102)). 
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expansion for ( |101| ) is that, in spite of the manifest linear 0{N) symmetry of ( |105| ), 
the nonhnearity of the "Goldstone" field n(x) is preserved on the quantum level as an 
identity on the correlation functions : 



Af\n'^P{n,dn) 



const (j\f P{n, dn) 



(106) 



where P(n, d-n) is arbitrary local polynomial of the fundamental fields and their deriva- 
tives, and A/'[. . .] indicates renormalized normal product of the corresponding composite 
fields. 

Armed with the above machinery, the first higher quantum conserved current (for 
s = 3 in the notations of (0), (0)) takes the following form : 



f w = ^^ 



+ aiAf 



n] 



e(2) 



+ a2]M\{d_n) 



a 



(107) 

where all coefficients ai,2,3 = 0{1/N) are expressed in terms of one-particle irreducible 
correlation functions and their derivatives in momentum space at zero external mo- 
menta. Their explicit form can be found order by order in 1/A^ from the renormalized 



1 / A^-diagram technique described above ||52 



Let us stress, that the higher quantum conserved currents ( p.07| ) and those for s 



5, 7, . . . do not have analogues in the classical conformally invariant theory |56 . 
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